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This paper begins with an examination of the revival structure and 
long-term evolution of Rydberg wave packets for hydrogen. We show 



that after the initial cycle of collapse and fractional/full revivals, which 
occurs on the time scale t rev , a new sequence of revivals begins. We 
find that the structure of the new revivals is different from that of the 
^ ■ fractional revivals. The new revivals are characterized by periodicities 

in the motion of the wave packet with periods that are fractions of 
the revival time scale t rev . These long-term periodicities result in the 
autocorrelation function at times greater than t Tev having a self-similar 
resemblance to its structure for times less than £ rev . The new sequence 
of revivals culminates with the formation of a single wave packet that 
more closely resembles the initial wave packet than does the full re- 
vival at time t rev , i-e., a superrevival forms. Explicit examples of the 
superrevival structure for both circular and radial wave packets are 
given. We then study wave packets in alkali-metal atoms, which are 
typically used in experiments. The behavior of these packets is affected 
by the presence of quantum defects that modify the hydrogenic revival 
time scales and periodicities. Their behavior can be treated analyti- 
cally using supersymmetry-based quantum-defect theory. We illustrate 



our results for alkali-metal atoms with explicit examples of the revival 
structure for radial wave packets in rubidium. 

To appear in the June 1995 issue of Physical Review A 
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I. INTRODUCTION 



A localized electron wave packet is formed when a short-pulsed laser excites a 
coherent superposition of Rydberg states |], @] • Wave packets of this type offer the 
opportunity to investigate the classical limit of the motion of electrons in Rydberg 
atoms. Initially, the motion of the wave packet follows the classical motion of a 
charged particle in a Coulomb field. The period of the motion is the classical period 
T c \ of a particle in a keplerian orbit. However, this motion persists only for a few 
cycles, whereupon quantum interference effects cause the wave packet first to collapse 
and then to undergo a sequence of revivals |]T], 0, |^, || . The revivals are characterized 
by the recombination of the collapsed wave packet into a form close to the original 
shape, which again oscillates with period T c \. The recombined wave packet is called 
a full revival, and it appears at a time t rev - For various times earlier than t vev , the 
wave packet gathers into a series of subsidiary wave packets called fractional revivals. 
The motion of these fractional revivals is periodic, with a period equal to a rational 
fraction of the classical keplerian period. 

This cycle of collapse, fractional revivals, and full revival is characteristic of sev- 
eral distinct types of Rydberg wave packets. Among these are radial, circular, and 
elliptical packets. A radial wave packet is localized in the radial coordinate but is in 
a definite eigenstate of the angular momentum. A circular wave packet is a super- 
position of fully aligned eigenstates with maximal values of the angular-momentum 
quantum numbers. It is localized in both angular and radial coordinates, and it fol- 
lows a circular trajectory. Lastly, an elliptical wave packet is localized in all three 
dimensions and travels along a classical keplerian orbit. 

Recently, there has been much experimental interest in the study of these Rydberg 
wave packets. Radial wave packets can be formed by the excitation of the atom from 
the ground state by a short-pulsed laser. The time evolution of the wave packet 
may then be studied using a pump-probe method of detection involving either time- 
delayed photoionization or phase modulation || ||, [|. The periodic motion of 
the radial wave packet, with the period equal to the classical orbital period, has 
been observed || and full and fractional revivals have been seen experimentally 
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Ill [12|]. Circular states have been produced as well, using either the adiabatic 



microwave transfer method |13| or crossed electric and magnetic fields ||14|| . These 
methods excite a range of states in a fixed n manifold, resulting in a wave packet 
that is stationary. A technique for the excitation of circular-orbital states to generate 
a circular wave packet has been proposed [EK]. If an additional weak electric field is 
present, an elliptical wave packet of arbitrary eccentricity could be produced. 

Several theoretical approaches have been used to study the properties of Ryd- 
berg wave packets. Some involve a description via different types of coherent states 
fill . Radial wave packets have been studied both numerically fl[] and perturbatively 
0, and recently a description as a type of squeezed state has been obtained |17|, [18| . 
They exhibit oscillations in the uncertainty product characteristic of a squeezed state. 
These states undergo full and fractional revivals with localized subsidiary waves dis- 
tributed along the radial direction. In contrast, the coherent-state approaches dealing 
with circular and elliptical states |19|, ^(], [21], ^2], [23], ^4], |2j| all involve superpositions 
of angular-momentum eigenstates and/or large angular- momentum quantum num- 
bers. The simplest such state is a circular wave packet consisting of a superposition 
of aligned eigenstates weighted by a gaussian function |2"5| . These have been studied 
numerically and exhibit both fractional and full revivals. The subsidiary waves that 
form in this case are at a fixed radius and are distributed in the azimuthal angle. 

Most work on the evolution of Rydberg wave packets has focused on the structure 
of the fractional revivals, which occur at times less than t rev . In this paper, we examine 
the evolution and revival structure of Rydberg wave packets for times beyond the 
revival time. We find that after a few revival cycles the wave packet ceases to reform 
at multiples of the revival time. Instead, a new sequence of collapses and revivals 
begins that is different from that of the usual fractional revivals. This sequence 
culminates with the formation of a wave packet that better resembles the original 
packet than does the full revival at time t iev . We refer to packets of this new type as 
super revivals. 

The appearance of the superrevivals is controlled by a new time scale t ST . We 
find that at certain times tf rac before t sx , the wave function can be written as a sum 
of macroscopically distinct subsidiary waves. Their motion is periodic with period 
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T{ ISLC given in terms of t iev and T c \. This periodicity is on a much longer time scale 
than that of the usual fractional revivals. For times greater than t rev , it produces in 
the autocorrelation function a self-similar resemblance to its appearance at times less 
than t rev . 

In addition to presenting the above qualitative features of the new revival struc- 
ture, this paper discusses a quantitative theory providing detailed expressions for the 
various time scales and periodicities that appear. The features of the theory are il- 
lustrated with explicit examples for both circular and radial wave packets. In these 
examples, we first consider wave packets having large values of the principle quantum 
number n, since these display detailed features of the superrevivals most strikingly. 
We then consider wave packets with smaller values of n for which the time scales for 
the superrevival structure fall within a range currently accessible in experiment. 

The above discussion applies to wave packets in hydrogen. In the latter part of 
this paper, we extend our analysis to obtain an analytical description for the case 
of Rydberg wave packets in alkali-metal atoms, typically examined in experiments. 
Examples are given to illustrate the resulting modifications of the behavior. The issue 
of distinguishing effects on the revival times and periodicities caused by the quantum 
defects from those caused by a detuning of the laser has been considered in ref. [EB 



In fact, not only are the effects different, but also the quantum-defect modifications 
to long-term revival times cannot be obtained by a direct scaling of the hydrogenic 
results. 

The structure of this paper is as follows. Sec. II presents the generic form we 
use for circular and radial wave packets and introduces the distinct times scales T c \, 
trev, and t sr determining wave-packet behavior over the time period of interest. In 
Sec. Ill, we examine the long-term behavior of Rydberg wave packets for hydrogen. 
We show that at certain times tf rac the wave function is well approximated as a sum 
over subsidiary waves with specified coefficients, and we present constraints on the 
allowed values of tf rac . We also demonstrate that at these times the motion of the 
subsidiary wave packets is periodic with a period T irac dependent on t Tev and T c \. The 
periodicity of the autocorrelation function is analyzed. Sec. IV provides illustrative 
examples of circular and radial wave packets in hydrogen both for large values of n 
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and for smaller values that are experimentally viable. 

In Sec. V, we extend our results to the case of wave packets in alkali-metal atoms, 
establishing the modifications in the revival time scales and periodicities that result 
from the inclusion of quantum defects. Explicit examples for the case of rubidium are 
provided in Sec. VI. Section VII discusses the results and considers some of the issues 
involved in experimental verification of our theory. Throughout the paper, technical 
details of the various proofs are relegated to appendices. 

II. TIME SCALES FOR WAVE PACKET EVOLUTION 

This section establishes some notation and conventions and provides a brief review 
of the evolution of Rydberg wave packets in hydrogen for times up to the revival time 
£ rev . In Sec. IIA, we present the wave function we use throughout the paper for 
the analysis. It is a generic form that permits simultaneous treatment of circular 
and radial wave packets. In Sec. IIB, we expand the component eigenenergies in a 
Taylor series about the wave-packet energy. This leads to the definition of distinct 
time scales governing the evolution of the wave function. Some background on the 
behavior of hydrogenic wave packets for times up to t rev is given in Sec. IIC, along 
with an example. 

A. Circular and Radial Wave Packets 

The time-dependent wave function for a hydrogenic wave packet may be expanded 
in terms of energy eigenstates as 



Here, E n = — l/2n 2 is the energy in atomic units, and i/j n i m (r) is a hydrogen eigenstate 
of energy and angular momentum. The weighting coefficients c n \ m = (VWi^l^O 5 ? 0)) 
depend on the initial wave function. 

Two types of wave packets used in experiments are considered in this paper. 
The first type is the circular wave packets, which consist of a sum of fully aligned 
eigenstates, i.e., ones with I = m = n — 1. In this case, the sum in Eq. ([I]) reduces 
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to a sum over n alone, and the weighting coefficients c n are independent of I and m. 
The second type is the radial wave packets, which consist of a sum of eigenstates with 
fixed angular-momentum quantum number. For example, all the eigenstates may be 
p states. Once again, the sum in Eq. (fj) depends only on n, and we may drop the / 
and m subscripts on c n . 

Adopting the generic notation that represents the eigenstates appropriate 

to a given type of wave packet, we may then rewrite Eq. (|T]) as 

oo 

t) = J2 c nVn(r) exp [-iE n t] . (2) 

n=l 

For the circular wave packets (p n {r) = ip n ,n-i,n-i(r) , while for the radial wave packets 
(Pn(r) = i^n,i,o(r). 

Both types of wave packet are excited by a short laser pulse. Since the laser can 
be tuned to excite coherently a superposition of states centered on a mean value of 
the principal quantum number n, in what follows we assume that the distribution is 
strongly centered around a value n. We may therefore approximate the square of the 
weighting coefficients as a gaussian function 

I |2 1 _ (n-n) 2 

\ C n \ — e 2 CT 2 (3) 

In the derivations that follow, we take n to be an integer. The results can read- 
ily be extended to incorporate noninteger n. The effects of noninteger n and the 



consequences of laser detuning are examined in ref. 26 



B. Time Scales 

If we expand the energy in a Taylor series around the centrally excited value n, 
we obtain 

E n ~ E n + E' n (n - n) + ±E%(n - nf + -E'^{n - nf + ■ ■ • , (4) 

2 o 

where each prime on E n denotes a derivative. The terms with derivatives define 
distinct time scales that depend on n. The first time scale, 

T cl = ^ = 2vr# , (5) 
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is the classical keplerian period. It controls the initial behavior of the packet. The 
second time scale, 

-27T _ 2n 

^rev — i rp„ ~ ~^~-*cI > (") 

is the revival time. It governs the appearance of the usual fractional and full revivals. 

The subject of this paper is the behavior of the packet on time scales greater than 
t iev . This behavior is dominantly controlled by a third time scale, 

2n 3n 

6 n ^ 

The scale t sr 3> t rev is a larger time scale we refer to as the superrevival time. Note that 
although t sr is large (typically about two orders of magnitude greater than t rev ) for the 
range of n of interest, it is still much smaller than the lifetime of the excited Rydberg 
atom. It is also smaller than the typical time scales for which microwave black- 
body radiation causes transitions [pTfl . It is therefore experimentally and theoretically 
reasonable to examine the behavior of the packet on this time scale. 

Keeping terms through third order, and defining the integer index k = n — n, we 
may rewrite Eq. (0) as 



k=—oo 



fkt k 2 t kH s 
-2ni + — 

V cl ^rev ^sr , 



(8) 



We assume that n is large so that the lower limit in the sum may be approximated 
by — oo. In this notation c k and (fkif) represent c n= k+n and <p n =k+n(r), respectively. 

C. Review of Behavior within the Revival Time Scale 

This subsection briefly summarizes the known behavior of the packet within the 
revival time scale and provides an example useful for comparison in our later analysis. 

The quantum evolution of the wave function (Q) up to times of order t rev is gov- 
erned by the first two terms in the exponential function in Eq. ([8]). If only the term 
linear in k were present, the wave packet would evolve like any localized packet for 
the harmonic oscillator, following the classical motion and oscillating indefinitely with 
period T c \. 
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The higher-order terms modify this behavior. The term quadratic in k governs 
the collapse and revival of the packet for times within t rev . An analysis of the role of 
this term is given in ref. ||, where it is shown that at the times t ~ ^t rev , where m 
and n are relatively prime integers, the wave packet may be rewritten approximately 
as an equally weighted sum of subsidiary wave packets. These are the usual fractional 
revivals. 

The behavior of the fractional revivals may be studied by examining the absolute 
square of the autocorrelation function Q, 

\A(t)\ 2 = \j2\cn\ 2 e- iEnt \ 2 . (9) 

n 

Since the wave function ip n (f) has been integrated over, we find the same autocor- 
relation function and hence the same revival structure for both circular and radial 
wave packets. At the fractional-revival times t ~ —t Iev , the autocorrelation function 
is periodic with fractional periods -T c \, where r denotes the number of subsidiary 
wave packets. 

An example useful for comparison with our analyses below is a wave packet in 
hydrogen centered about n = 319 with a = 2.5. Figure 1 shows the square of the 
autocorrelation function for times up to and just beyond t rev ~ 1.05 //sec. Fractional 
revivals are prominent at t ~ |i rev — 0.13 //sec, t ~ |t rev — 0.18 /tsec, t ~ ^t vev ~ 0.26 
/tsec, t ~ |t rev — 0.53 //sec, and at multiples of these values. The periodicity of the 
autocorrelation function at the one-half and full revivals is T c \ ~ 4.9 nsec. This can 
be seen in the figure. The periodicities at the earlier fractional revivals, for t ~ |t rev , 
gi rev , j^rev, respectively, are jT c i, |T c i, \T C \. As can be observed in the figure, the 
peaks in the autocorrelation function are not as pronounced for times approaching 
trev as they are for earlier times. This is because higher-order terms in the energy 
expansion ([|) act to distort the revivals. 

III. WAVE PACKETS IN HYDROGEN 

In this section, we show that at certain times tf rac it is possible to expand the 
third-order wave function ty(r,t) of Eq. (§) as a series of subsidiary wave functions. 
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Following ref. ||, the idea is to express \l/(r, t) as a sum of wave functions ip c \ with 
matching periodicities and with shape similar to that of the initial wave function 
\I/(r, 0). The form of the coefficients in the expansion of ^f(r, t) is obtained and used 
to show that at times tf rac the expansion correctly reproduces the wave function. At 
particular values of if rac , superrevivals are shown to form that more closely resembles 
the initial wave packet than does the full revival. Lastly, we show that the wave 
function and autocorrelation function are periodic with a given period Tf rac . 

In what follows, we find it is useful to write the times if rac at which the subsidiary 
wave packets form as a linear combination of rational fractions of t ST and t rev . We 
therefore define 

V m / 1n \ 

£frac — ^sr ^rev > (1") 

q n 

where p and q are relatively prime integers, as are m and n. Constraints on the 
possible values of these integers are derived below. In particular, we find that m is 
nonzero whenever n cannot be evenly divided by four. We therefore write n as 

n = Ar] + \ , (11) 

where 77 and A are both integers, and A = 0, 1, 2, or 3. 

The subsidiary wave functions we use in the expansion of the full wave function 
are defined in terms of a function tp c \ depending only on the first-derivative term in 
the energy expansion, 

oo 

ip cl (f,t)= c k Lp k (f)exp 

k=— oo 

We write it in terms of t vev instead of T c \ because the former is more useful when 
considering times greater than t Tev . 

The subscript on ip c \ is a reminder that this wave function follows the classical 
motion at all times. Note also that at time zero, ijj c \(r,0) = \I/(r, 0), i.e., tp c \(f,0) 
matches the initial wave function exactly. These features of the function i)j c \(r,t) 
make it a suitable candidate for use in generating an expansion of the full wave 
function *&(r, t) as a sum of macroscopically distinct subsidiary wave functions. 

Appendix A shows that at the times tf rac the wave function $?(r, t) can be written 
as a sum of subsidiary waves having the form of ipdif, t), but at times that are shifted 



2n kt 
-2ni 

3 trev 



(12) 
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by a fraction of the period T c \. The shifted functions have the form ip c i(r,t + ^fT c \), 
where / and a are certain specified integers and s = 0,l, 2, — 1. Since these 

functions form a set with the same periodicity as ty(r,t) at the times tf rac , we can 
use them as a basis for expanding \l/(r, t). Moreover, the shapes of the shifted ip c \ 
resemble the initial wave function at different points in its cycle, after it has been 
time-translated by a fraction of the classical period T c \. An expansion of this type 
provides a natural formalism within which to examine the occurrence of long-term 
revivals. 

We therefore write an expansion in subsidiary wave functions at times t ~ tf rac in 
the form 

l ~ x sa 
^(r,t) = J2bsMr,t + —T cl ) . (13) 

s=0 1 

The coefficients b s are complex valued and are given as 

exp [2m9 k/ ] , (14) 



1-1 as 



'■- - tE ex P 



^ fc'=o L I 
where 8^ is defined in Eq. ( |3~4"D of Appendix A. 

Appendix A also derives the constraints on the integers p, q, m, and n in the 
definition fllPf) of £f rac , and shows that p — 1, q is a multiple of three, and m/n obeys 
Eq. (|33|). With these constraints, the expansion of ^(r, t) in terms of subsidiary 
waves ip c \(r,t) holds at times £f rac given by Eq. (|T0|). Since t sr 3> t rev , we conclude 
that interesting behavior of the wave packet can be expected near times that are 
simple fractions of |t sr . This somewhat counterintuitive result is strongly supported 
by the examples presented in Sec. IV. 

Whereas the usual fractional revivals consist of r subsidiary packets weighted 
equally by factors of 1/r, and hence have fractional period T = T c \/r, we find the 
superrevivals have different behavior. An interesting feature of the b s coefficients (14) 
at the times tf rac is that they may have different moduli |6 S | 2 , and hence distinct 
subsidiary wave packets may not all have the same weight. Explicit examples of this 
are given in Sec. IV. In certain cases, all the b s coefficients except one vanish. The 
remaining one then has modulus one, corresponding to the formation of one wave 
packet. It turns out that at the times tf rac , even if the subsidiary wave packets are 
unequally weighted, the wave packet and the autocorrelation function are periodic. 
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In what follows, it is useful to introduce terminology distinguishing the different 
types of superrevival that can appear. We refer to the set of subsidiary packets at 
the times tf rac as a fractional superrevival. However, occasionally only a single packet 
appears, resembling the initial wave packet more closely than the usual full revival 
does at the time t Tev . This is because contributions are incorporated from higher- 
order corrections in the expansion of the time- dependent phase. We call a single 
wave packet of this type a full superrevival. 

After the formation of a fractional superrevival consisting of several subsidiary 
wave packets, it often happens that the subsidiary packets quickly evolve into a 
configuration where one of them is much larger than the others. The dominant wave 
packet in this case can again resemble the initial wave packet more closely than does 
the full revival at time t rev . We refer to a configuration of this type, consisting of 
primarily one large wave packet, as a partial superrevival. 

Appendix B proves that the wave packet \^(f, t)\ 2 as well as the absolute square 
of the autocorrelation function A(t {rSiC ) = (^(r, 0)|\l/(r, £f rac )) are periodic at times 
t ~ tfrac- The period T frac is given as 

Tfrac ^rev T c \ , (15) 

q V 

where the integers u and v satisfy Eq. (|60|) . In addition, the value of p in t frac must 
equal 1 for this periodic behavior to occur. 

This periodicity is different from that of the usual fractional revivals, for which 
the periods are fractions of T c \. Instead, the periods of the fractional superrevivals are 
fractions of t rev combined with small additional shifts that depend on T c \. The time 
scale is considerably greater for the fractional superrevivals than for the fractional 
revivals. In addition, this periodicity causes the autocorrelation function for times 
much greater than t vev to have a self-similar resemblance to its form for times less 
than t rev . We will illustate this behavior explicitly in the next section. 

IV. EXAMPLES FOR HYDROGEN 

This section presents several examples illustrating in detail the various aspects of 
the long-term evolution of Rydberg wave packets in hydrogen. Results numerically 
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generated from analytical expressions are compared with the theoretical predictions 
given above. 

In these examples, we assume the laser is tuned to excite coherently a distribution 
strongly centered at either fi = 319or?7, = 36. The former is large and ensures that 
interesting features of the the superrevival structure are noticeable. Since 319 is not 
divisible by four, the quantity A in Eq. ([□]) is nonzero, which in turn makes the 
ratio m/n nonzero in the definition of £f rac . The full analysis of Sec. Ill is therefore 
required. An example with n = 320, which is divisible by four, can be found in ref. 

The other choice, n = 36, is motivated by experimental considerations. For this 
case, the associated time delay in a pump-probe experiment falls within a range that 
is currently accessible. Our analysis shows that the corresponding full superrevival oc- 
curs after approximately 776 psec, which is within the range of time delays considered 
in refs. |TO|, [TT|. 

Sec. IVA examines the autocorrelation functions for the example with n = 319. 
Since the autocorrelation functions for corresponding circular and radial wave packets 
are identical, this provides a relatively broad perspective on the long-term behavior. 
In Sec. IVB, we examine detailed features of circular wave packets with n = 319 at 
the times tf rac . Hydrogenic radial wave packets with the experimentally accessible 
value n = 36 are considered in Sec. IVC. 

A. Autocorrelation Functions 

This subsection discusses the autocorrelation functions for our two choices of n. 
The absolute squares of the autocorrelation functions are computed directly using the 
analytical expression in Eq. (^) and the definition of the c n coefficients given in Eq. 
(0). These results are then compared with the preceding theoretical analysis. In what 
follows, we restrict ourselves to times of order |t sr , corresponding to a minimum value 
of q equal to six. An accurate treatment on longer time scales requires fourth-order 
and higher terms in the Taylor-series expansion of the energy. 

The autocorrelation function is primarily useful in the evaluation of theoretical 
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predictions because it is sensitive to the formation of single packets. Thus, in the 
following examples both full and partial superrevivals appear as distinct peaks. How- 
ever, fractional ones are less apparent. For this reason, we defer to Sees. IVB and IVC 
comparison of the results of numerical computations with any details of our theory 
involving, for instance, the explicit values of the coefficients b s . 

As an example, consider the autocorrelation function for a Rydberg wave packet 
in hydrogen with n = 319 and o = 2.5. In this case, rj = 79, A = 3, and t ST ~ 251 
/isec. The absolute square of the autocorrelation function is shown in Fig. 2. Figure 
2a shows the behavior of the wave packet for the first 15 /isec. The first 1.2 /xsec 
contains the cycle of fractional and full revivals shown in Fig. 1. This cycle collapses 
after approximately 5 /isec, whereupon new structure emerges. Figures 2b and 2c 
show |A(£)| 2 for times up to and just beyond |t sr . 

Our prediction for q = 6 is that a single wave packet should form at t « tf rac ~ 41.4 
/xsec. Only one wave packet appears because only one h s coefficient is nonzero (see Sec. 
VIB). Equation (^) predicts the period of this packet as Tf rac ~ 0.52 //sec. In Fig. 2c, 
we observe large peaks in |^4(t)| 2 near 42 /zsec, with a periodicity approximately equal 
to Tf rac . These peaks are larger than the full-revival peak at t ~ t rcv ~ 1.05 /zsec. This 
is because at t ~ 41.4 /xsec the higher-order corrections to the time-dependent phase 
are contributing coherently to the wave function. Note also that between the large 
peaks at t ~ 41.4 //sec there are smaller peaks with periodicities that are fractions 
of Tf rac . These correspond to the formation of fractional revivals during the times 
between the periodic appearances of the full superrevival. 

For q = 9, the theory gives tf rac ~ 27.6 /isec and Tf rac ~ 0.35 [isec. Fig. 2b shows 
that the autocorrelation function exhibits peaks near this time, but with approxi- 
mately half the period Tf rac . This occurs because more than one subsidiary wave 
packet is forming, and they are unequally weighted. In this case, the fractional su- 
perrevival evolves into a partial superrevival at two separate times during each cycle. 

The q = 12 peaks at tf rac ~ 20.7 //sec with Tf rac ~ 0.26 /zsec are more pronounced 
than the q = 9 peaks. This is because the corresponding subsidiary wave packets 
evolve into a structure with a single dominant wave packet, thereby creating large 
peaks in the autocorrelation function. The periodicity of the peaks agrees with the 
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predicted value of Tf rac . 

No prominent peaks appear in Fig. 2b for q = 15 at the times tf rac ~ 16.6 /isec. 
As we show in the next subsection, in this case there are four unequally weighted 
subsidiary wave packets. These evidently do not evolve into a structure with one 
dominant wave packet. 

Since q = 18 is divisible by 9, the number of terms in the expansion of the wave 
function in this case is I — | = 6. Of these six terms, three vanish. The remaining 
terms are unequally weighted with one of them substantially larger than the rest. 
This leads to the formation of a single packet at time tf rac ~ 13.8 /isec, with a period 
Tf rac ~ 0.17 /isec. Figure 2a shows peaks in the autocorrelation function at this time 
with a periodicity that agrees with Tf rac . 

Peaks also appear in the autocorrelation function near t ~ 7 /isec. For q = 36, the 
theory gives tf rac ~ 6.91 /isec and T frac ~ 0.09 /isec, which agrees with the observed 
periodicity. In this case, however, we find that there are 12 nonvanishing b s coefficients 
that are unequal in magnitude. It is unclear why these subsidiary waves evolve into 
a partial superrevival. 

B. Circular Wave Packets 

This subsection examines more closely our theoretical description of circular Ry- 
dberg wave packets in hydrogen with n = 319 and a = 2.5. The time-dependent 
wave function is given analytically in Eq. (0) with the c n coefficients specified by Eq. 
(0) and with y3 n (r) = ip n , n -i tn -i{r). Plots of the wave packet at different times can 
be obtained numerically by taking advantage of the peaking of the c n coefficients 
around the value n. Figure 3 shows the azimuthal dependence at various times of a 
cross-sectional slice of the wave packet taken in the plane of the orbit and at fixed 
radius equal to the expectation value (r) = \fi (2n + 1). 

The initial wave packet at time t = is presented in Figure 3a. It is localized 
about the azimuthal angle = 0. Figure 3b shows the first full revival, at time t rev . 
At this point, the packet has collapsed, passed through the sequence of fractional 
revivals, and reformed. It is primarily a single packet shifted in phase relative to the 
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initial one, but some smaller subsidiary packets remain visible. 

The remainder of the graphs in Fig. 3 show the wave packet at times relevant for 
comparison with our theoretical analysis. In what follows, we give the b s coefficients 
for the expansion (|T3|) of the wave packet at the times £f rac for each of the cases q — 6, 
9, 12, 15, and 18. We also compare the packets at the times tf rac and tf rac + Tf rac to 
test the predicted periodicity. For n = 319, we have rj = 79 and A = 3. The value of 
I is given in Eq. and depends on q. The integer a that appears in the expansion 
of \P(r, t) in Eq. (|I~3"D is equal to 8r]/N where N is the product of all the factors of 
(877) that are also factors of I. Note that we set the integer a in Eq. (fiTf ) equal to 1. 
In fact, its value is irrelevant: if it changes, the b s coefficients become permuted but 
leave unaffected the total sum of the subsidiary waves. 

For q = 6, we obtain / = 6, N = 8, a = 79, tf rac = |t sr — |t rev , and T frac = 
\trev — \T C \. Using Eq. (|14D , we find that 64 = 1, while all the other b s coefficients 
vanish. The wave function at tf rac can therefore be written as a single subsidiary 
packet. Since ipci is periodic with period T c \, it is sufficient to evaluate the phase 
of ?/>ci(£frac + x T c \) for s = 4 modulo T c \. Suppressing the r dependence, we find 
^(£frac) ~ ^ci(O). Figure 3c shows the wave packet at the time £f rac . It does indeed 
consist of a single wave packet at the initial point in the orbital cycle. Comparing 
Fig. 3c to Fig. 3b, which represents the wave packet at the full revival, we see that 
the full superrevival resembles the initial wave in Fig. la more closely than does the 
full revival at t — t vcv . Moreover, the tail of subsidiary waves visible in Fig. 3b is 
absent in Fig. 3c. The wave packet at the time tf rac + Tf rac is plotted in Fig. 3d. This 
figure resembles Fig. 3c, which verifies the predicted periodicity of the wave packet. 

With q = 9, we obtain I — 3, N — 1, a — 632, tf rac = |t sr — ^t rev , and T frac = 
l^rcv — f^ci- I n this case, there are three nonzero terms in the expansion of the wave 
function. We find 

#(*frac) » b iJ cl (^) + M> d (~) + &2^cl(^) , (16) 

where the argument of ipci is written in units of T c \, and where |&o| — 0.45, ~ 0.85, 
and I&2I — 0.29. This expansion predicts that the wave packet consists of three 
subsidiary packets. Since t/'cI passes through the range to 2n in at a constant rate, 
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at a specified fractional interval in the cycle its angular position is at the corresponding 
fractional part of 2ir. We therefore expect a large wave packet at ~ |(27r) , a medium 
one at ~ §(27r), and a small one at ~ §(27r). This prediction agrees with the 
plot in Fig. 3e. The wave packet at the time tf rac + Tf rac is shown in Fig. 3f. The 
periodicity Tf rac of the wave packet is evident. 

Since the wave function in this case has three distinct subsidiary wave packets, a 
large peak in the autocorrelation function is not to be expected. This is compatible 
with the structure of Fig. 2b at times near t frac ~ 27.6 /jsec. The appearance of peaks 
of equal height with a periodicity |Tf rac indicates that the fractional superrevival is 
evolving into partial superrevivals twice each cycle, which is possible because the 
subsidiary wave packets are unequally weighted. 

For q = 12, we find I = 12, N = 8, a = 79, t frac = ^t sr - ^t rev , and T frac = 
\t rev — jjTd. The expansion of the wave packet yields 

*(*frac) « Mcl(Jj) + M>d(^) + M>d( J) + 6io^d(0) (17) 

with \b\\ = 1 64 1 = I b 7 1 = |6 10 | — \. The remaining b s coefficients vanish. The packet 
in this case consists of four equally weighted subsidiary packets evenly distributed in 
0. Figure 3g confirms this prediction. However, there is some distortion, arising from 
higher-order terms neglected in our analysis. In Fig. 3h we plot the wave packet at 
time tf rac + Tfrac- The periodicity is again evident. 

For q = 15, we have I = 15, N = 1, a = 632, t frac = ^t sr - ^t rev , and T frac = 
|t rev — y|T c1 . The expansion gives 

* (tfrac) ~ 62^(5) + hM^) + bsM^) + 6i4^d(0) , (18) 

where the moduli of the complex coefficients are I&2I = |&s| — 0.45, |6g| ~ 0.72, and 
|6 14 | ~ 0.28. For these values, we expect two equally weighted subsidiary wave packets 
at §(27t) and |(27r), a large subsidiary packet at |(27r), and a small one at the origin. 
Figure 3i displays the packet at tf rac . It agrees with the predictions. Figure 3j shows 
the wave packet at the time tf rac + Tf rac , which again exhibits the expected periodicity. 

Finally, in Fig. 3k, we show the wave packet corresponding to q — 18. Here, we 
obtain I = 6, iV = 8, a = 79, tf rac = j^t ST — |t rcv , and Tf rac = |t rev — j|T cl . There are 
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three nonzero coefficients b s , so 

14 7 

#(tfrac) ~ boM-) + hMg) + M^i(-) , (19) 

with |&o| — 0.84, 1 2>2 j ~ 0.29, and |£>4| ~ 0.45. A large subsidiary wave packet is visible 
near h(27t), an intermediate-sized one is near §(27r), and a small one is near |(27r), in 
agreement with the theoretical result. Figure 31 shows the wave packet at tf rac + T frac , 
confirming that the periodicity is correctly given by Tf rac . 

For q greater than 18, more than four subsidiary wave packets appear in the 
expansion of the wave function at the times tf rac . There are therefore fewer pronounced 
peaks in the autocorrelation function, and prominent partial superrevivals are less 
likely to form. 

C. Radial Wave Packets 

The previous sections have established the validity of our theory for the various 
types of superrevival. This subsection treats an example closely related to a feasible 
experimental situation: a p-state radial Rydberg wave packet for hydrogen with n = 
36 and a = 1.5. 

Figure 4a shows the plot of the radial wave packet at time t = 0. The wave packet 
starts its cycle at the inner apsidal point and then moves radially to the outer apsidal 
turning point. Initially, the wave packet is highly oscillatory. As it moves, the packet 
changes shape in a way characteristic of a squeezed state |T7|]. In Fig. 4b, we show 
the wave packet at the outer apsidal point at the time \T C \. Here, it is localized in 
the radial coordinate and is at its point of minimum uncertainty. 

At the full revival, t = t rev , the wave packet has reformed after undergoing its 
cycle of collapse and fractional revivals and is located at its initial position near the 
inner apsidal point. For ease of comparison of the initial and full- revival packets, we 
show in Fig. 4c the full- revival packet displaced forward in time by \T C \. This brings 
it to the outer apsidal point where it is the most localized. The initial and full-revival 
packets are similar, although the latter has some distortion due to the higher-order 
corrections. 
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As before, the full superrevival is predicted to occur for q = 6. This gives I = 6, 
N = 72, a = 1, tf rac = ht ST , and T frac = ~t rev — ~T c i. Note that since n = 36 
is divisible by four, the fraction m/n in tf rac vanishes. The expansion of the wave 
function is ^(tf rac ) ~ ^4"0ci(|)j with 64 = 1 and all the other b s coefficients vanishing. 
In this case, ipci is a radial wave function, which in one complete cycle moves from the 
inner turning point to the outer one and back. At tf rac , the wave function therefore 
represents a single packet two-thirds through the classical cycle. To compare this 
packet to the full revival, we view it displaced backward by ~T c i so it is at the outer 
turning point. Figure 4d shows the wave packet at this time. The full superrevival is 
visibly better localized than the full revival in Fig. 4c, and it has less distortion. In 
Fig. 4e, the wave packet is displayed one period Tf rac later than the time of Fig. 4d. 
Comparison shows that the packet has the predicted period T frac . 

V. RYDBERG WAVE PACKETS IN ALKALI-METAL ATOMS 

The previous sections of this paper have provided a satisfactory theory of the 
long-term behavior of Rydberg wave packets in hydrogen. However, the wave packets 
used to date in experiments to exhibit full and fractional revivals || || |H], [11], H3 
have been produced in alkali-metal atoms. In the remainder of this paper, we show 
how to extend our previous theoretical analysis to these cases. 

The Rydberg series for an alkali-metal atom is given by the energies 

where n* = n — 5(1), and 6(1) is an asymptotic quantum defect for an alkali- Rydberg 
atom. In what follows, we consider Rydberg wave packets in alkali-metal atoms that 
can be represented as a superposition of states strongly peaked around a central value 
n* — n — 6(1) of the shifted principal quantum number. Following the simplification 
introduced at the end of Sec. IIA, we assume n is an integer. The extension of the 
results below to noninteger n, along with the issue of disentangling laser detuning 



from features induced by quantum defects, has been considered in ref. [26]. 
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We proceed by expanding the Rydberg energies as a Taylor series in n*, 

E n * - En* + K* (n - n) + \e% (n - nf + ijgft (n - nf + • • • . (21) 

The derivatives are taken with respect to n*, and we have used the relation (n* — n*) = 
(n — n). The first three derivative terms define the time scales 

T* = 27m* 3 (22) 
2n* 

C = — T cl (23) 

and 

3n* 

^sr = ^^^rev J (24) 

which depend on the quantum defects. 

Since the quantum defects 5(1) depend on /, we may no longer treat circular and 
radial wave packets generically. Circular wave packets have large values of I and 
so vanishing values of 5(1). A hydrogenic description therefore should be a good 
approximation for the analysis of the behavior of circular wave packets for alkali- 
metal atoms. In contrast, radial wave packets have small values of /. This means 
their behavior is modified by the presence of quantum defects, with time scales and 
periodicities changed relative to those in hydrogen. In the following, we limit ourselves 
to considering radial wave packets, for which the quantum defects are important. 

To incorporate non- hydrogenic features in our treatment, we use supersymmetry- 
based quantum-defect theory (SQDT), which has analytical wave functions with exact 



eigenvalues reproducing the Rydberg series for alkali- metal atoms | 29fl . This model 
describes the behavior of the valence electron in an alkali-metal atom as that of a 
single particle in an effective central potential. The effective potential is found by 
acting on the radial Coulomb potential for hydrogen with a supersymmetry trans- 
formation and then adding specific supersymmetry-breaking terms that incorporate 
electron-electron interactions and reproduce the Rydberg series while leaving the 
eigenfunctions analytical. A modified angular quantum number I* = I — 5(1) + 1(1) is 
introduced, where 1(1) is an integer playing the role of the supersymmetric shift. The 
resulting SQDT hamiltonian is equivalent to one obtained by replacing n, I, E n in the 
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radial equation by n*, /*, E n *, respectively. The exact three-dimensional SQDT wave 
functions are Yi m (9,(f>)R n *i*(r). These have the eigenenergies E n * given in Eq. ( pT|) 
and form a complete and orthonormal set. More details on this model and references 



to the recent literature may be found in the review article [30 



Let us denote by </9* n (r) the eigenstates of an alkali-metal atom that are appropri- 
ate for the description of a radial wave packet. The analogue of Eq. ([I]) for the wave 
function of the packet is then given by 

oo 

*(^> *) = 2 c n *y?* n (r) exp \-iE n *t] , (25) 

71=1 

where the c n * are complex weights determined by the initial wave function. One role 
of SQDT in the analysis that follows is to provide explicit analytical expressions for 
the eigenstates in Eq. (p5|) , i.e., we take (p* n (r) = Yi m (0,<j))R n *i*(r). 
Using the Taylor series (|2T|), the wave function becomes 



^{r,t) = c fc v3*fc(f)exp 

k=— oo 



fkt kh k 3 t" 

-2ni - + — 

\ ± c\ b rev L sv , 



(26) 



We have kept the first three terms in the Taylor expansion so that times of order t* r 
can be considered. The sum ranges over the integer values of k = (n* — n*) = (n — n). 
We assume that n* is large so that the lower limit in the sum may be written as — oo. 
Since k is integer valued, we can take the coefficients Ck as the gaussian functions 
defined for hydrogen in Eq. @. 

We parametrize the times at which subsidiary wave packets form as 

t* — -t* — — t* (97) 

''frac _ g sr ^ ''rev > ' J 

where p and q are relatively prime integers, as are m and n. A priori, the integers p, 
q, m, n are unconstrained, and in particular they are not related to the corresponding 
quantities for a hydrogenic radial wave packet. 

Following the hydrogenic notation as closely as possible, we express n* as 

n* = + \ - - , (28) 
v 

where r], A, fi, v are all integer valued, A = 0, 1, 2, or 3, and \ijv is an irreducible 
fraction less than one. In the most general case, [ijv would represent the fractional 
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part of a combination of the quantum defect and a laser detuning. In this paper, we 
assume the laser is on resonance, so \ijv is the fractional part of the quantum defect. 



The effects of an additional laser detuning are described in ref. |28]. Note that writing 
the quantum defect as an integer part minus a rational fraction of unity can always 
be done for experimentally determined quantum defects, which of necessity have only 
a finite accuracy. 

We also define a wave function with time dependence involving only the linear 
term in k. In terms of t*. it is 



rev' 

2n* kt 



-2ni- 



''rev 



(29) 



ip c \(r,t)= c fc y?* fc (r)exp 

k=— oo 

Note that at time zero, \P(r, 0) = ^ c i(f, 0). 

In Appendic C, we show that the phases induced by the shifts At in T/> c i(r, t) have 
the same period as the higher-order contributions to the time-dependent phase in 
\l/(r, t) at the times t* h&c . It is therefore plausible to use the set ^ci(^£f rac + x^ci) 
with s = 0, 1, — 1 as a basis for an expansion of ^(f, tf rac ). We write 

i=J sa 

W, <£J = E bsMr, tfrac + -T^cl) • (30) 
s=0 f 



It is shown in Appendix C that the b s coefficients have the same form as in Eq. ( |14[ 
for hydrogen. The proof that this expansion is valid then follows as before, using Eq 



Equation (|30| ) implies that the expansion in terms of subsidiary packets of a wave 
packet in an alkali-metal atom is similar to a corresponding one in hydrogen. In 
particular, the number of subsidiary wave packets that form at times tf rac and their 
relative proportions are the same as those in a corresponding wave packet in hydrogen 
at times tf rac . However, the time scales t* r , t* ev , and T c *j controlling the behavior are 
different from those in hydrogen, as are the times tf rac at which the expansion in 
subsidiary packets is valid. 

In Appendix C, we determine the allowed values of tf rac . We find q must again be 
restricted to multiples of 3. The ratio m/n is given in Eq. fl67|). 

The proof that the wave packet and absolute square of the autocorrelation function 
are periodic for times near t£ rac , with p — 1, is also outlined in Appendix C. The 
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allowed values of the period T f * ac are found to be 

rp* +* rp* 

J frac — g rev y cl ' ^ ' 

where the integers it and v are given in Eq. (|75|). 

Comparison ofu/v for alkali-metal atoms (Eq. ( |75| ) of Appendix C) with the cor- 
responding definition for hydrogen (Eq. ( |6~0l) of Appendix B), reveals the appearance 
of an additional shift that depends on the fractional part fi/u of n*. This means 
that the quantum defects cause an extra shift in the period T f * ac , in addition to the 
rescaling of t rev and T c \. Note that if \x = 0, so that n* is an integer, then the analysis 
reduces to the hydrogenic case. 

To summarize, we have shown that a radial Rydberg wave packet in an alkali- 
metal atom forms superrevivals at times £f rac given by Eq. (|27|), where p — 1, q is a 
multiple of three, and the ratio m/n is determined. At these times, the wave function 
gathers into a finite number of subsidiary wave packets that move periodically with 
period T f * ac given by Eq. (|31|) with a specified ratio u/v. The expressions for tf rac and 
T f * ac differ from those for a corresponding wave packet in hydrogen by rescalings and 
shifts that depend on the quantum defect. 

VI. EXAMPLES FOR RUBIDIUM 

To illustrate in an experimentally viable scenario the long-term behavior of Ry- 
dberg wave packets when quantum defects are present, we consider in this section 
an example of a radial wave packet in rubidium with n = 36. For definiteness, we 
take p-state angular distributions, corresponding to the wave packets produced by a 
single short laser pulse. The example illustrates that full and fractional superrevivals 
may be observed for values of n corresponding to time delays currently accessible in 
experiments. Some other experimental issues are considered in the next section. 

In addition to providing an additional check on the predictions of our theory, the 
long-term revival structure arising in this example can be compared to the analogous 
situation for hydrogen presented in Sec. IVC For the autocorrelation function the 
behavior is generated from Eq. @ with E n replaced with E n *. For the radial wave 
packet itself, it is obtained from the analytical expression (|25|). We use a gaussian 
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distribution of width a = 1.5 for the weighting coefficients c n , which makes it a good 
approximation to truncate the sum over the SQDT eigenstates after a finite number 
of terms. 

The quantum defect for p states of rubidium is 5(1) ~ 2.65. This gives n* ~ 33.35, 
7] = 8, A = 2, and \xjv = 13/20. The time scales determining the behavior of this 
packet are therefore T c * ~ 5.63 psec, t* ev ~ 0.12 nsec, and t* T ~ 3.13 nsec. 

Figure 5 shows the absolute square of the autocorrelation function for this wave 
packet on time scales of order |t* rac . The full revival peak is visible near £* ev . For the 
relatively small value of n considered here, the revival structure is not as prominent 
as for larger values. Full and fractional superrevivals should nonetheless occur, albeit 
with less visible features. For q = 6, our theoretical analysis predicts the occurrence 
of a full superrevival at i* rac ~ 0.50 nsec with periodicity T f * ac ~ 0.06 nsec. Figure 5 
contains peaks near £* rac in agreement with these values. The theory also predicts the 
appearance of fractional superrevivals for q = 9 and 12 at the times tf rac ~ 0.33 nsec 
and 0.25 nsec, respectively. The corresponding predicted periodicities are T f * ac ~ 0.04 
nsec and 0.03 nsec. Peaks with appropriate periodicities may be observed near these 
times tf rac in Fig. 5. 

Consider next the radial wave packet itself for this example. Fig. 6 displays the 
packet at different times. The results should also be compared with the figures for 
the corresponding hydrogenic packet shown in Fig. 4. 

The initial packet at time t = is displayed in Fig. 6a. It is highly oscillatory 
and located near the inner apsidal point. Halfway through the classical orbit, at 
time t = IT*!, the wave function is radially localized and at its point of minimum 
uncertainty. See Fig. 6b. The wave packet near the full revival time t Tev is shown in 
Fig. 6c, after a time displacement by £T£. This displacement brings the packet to 
the outer apsidal point for ease of comparison. 

The theoretical analysis suggests that a full superrevival occurs at £* rac = |t* r — 
j^t* cv with period T f * ac = |t* ev - §jT c *. The definition ([14]) of the b s coefficients gives 
bi = 1, while the other b s vanish. Evaluating the expansion (p0[) yields \l/(tf rac ) ~ 
V'ci(^)) where the argument of ip c \ is written in units of T* v For purposes of compar- 
ison at the outer apsidal point, we displace the packet by an additional amount |T c *j. 
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The result is shown in Fig. 6d. A single wave packet appears, with less distortion 
than the full revival in Fig. 6c. 

In Fig. 6e, the wave packet is displayed a time T f * ac later than in Fig. 6d. The 
wave packet is again at the outer apsidal point. The motion has period in agreement 

with T f ; ac . 

VII. DISCUSSION 

In this paper, the behavior of Rydberg wave packets has been discussed on time 
scales large compared to the full-revival time t rev . We have demonstrated the ap- 
pearance of various types of superrevival and have presented a theoretical analysis 
that correctly predicts the times of formation, the associated periodicities, and their 
structure. Our analysis covers both hydrogenic Rydberg wave packets and the exper- 
imentally preferred Rydberg packets in alkali-metal atoms. 

The fractional superrevivals that we present have long-term periodicities different 
from those of the usual fractional revivals. These periodicities could be observed for 
values of n corresponding to time delays in a pump-probe experiment of under one 
nanosecond. 

When several packets form, we have shown that they often evolve quickly into a 
state where one packet is much larger than the others. In this case, large peaks appear 
in the autocorrelation function. All these features have been confirmed in numerically 
computed examples both for the autocorrelation function and for the wave packets 
themselves. 

The wave-packet behavior on the time scale t sr is self-similar in some ways to the 
fractional-revival structure appearing on the time scale t rev - There are significant and 
experimentally observable differences, however. For example, the times tf rac are more 
restricted than the times at which the usual fractional revivals occur, since only times 
approximately given by certain irreducible fractions of |t sr are permitted. Moreover, 
the allowed periodicities Tf rac are dominated by multiples of 3t rcv . The origin of the 
frequent appearance of multiples of three and one-third in our theory can be traced to 
the appearance of the third-order derivative of the energy in the original wave-packet 
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expansion. 

Our generalization of the hydrogenic analysis to include quantum defects uses a 
supersymmetry-based quantum-defect theory that provides a complete and orthonor- 
mal set of analytical wave functions with exact Rydberg eigenenergies. We have 
thereby proved that an expansion in terms of subsidiary packets is valid for a radial 
Rydberg packet in an alkali- metal atom. The controlling time scales T c \, t rev , and t sr 
are modified relative to the hydrogenic case by quantum defects, as are the periodic- 
ities Tf rac of the motion and the times £f rac at which the subsidiary-packet expansion 
is valid. Since the quantum defects for p states of the heavier alkali-metal atoms can 
be relatively large, these modifications can be significant and must be included in any 
accurate description of the evolution of the wave packet. 

Relatively little prior literature exists on the long-term behavior of Rydberg wave 
packets. A long-term revival at t — n 3 T c i has been presented in ref. Jl|]. The time- 
dependent autocorrelation function (Q) is an example of an almost-periodic function 
if the sum is truncated to a finite number of terms. Discussions of such functions exist 
in a more general context (see, for example, ref. However, the wave function 

expansion (EST) has spatial dependence that is critical in the superrevival structure, 



involving the formation of the subsidiary wave packets. Ref. discusses a class 
of long-term revivals in hydrogen occurring when t and n are such that the energy 
expansion yields a phase that is an integer multiple of 2iri. In contrast, the analysis of 
the full/fractional superrevivals presented here, which holds for arbitrary values of n 
in any alkali- metal atom, is based on finding the times tf rac where the phase in ^/(f*, t) 
matches the phases of the shifted waves ip c \ in the expansion. We find that at tf rac for 
q = 6, corresponding to the full superrevival, the phase of the time-dependent terms 
is not necessarily an integer multiple of 2iri. The point is that the superrevival times 
we find correspond to the formation of subsidiary wave packets. Unlike the usual 
fractional revivals, these are not equally weighted. The counterintuitive appearance 
of the q = 6 superrevival occurs because all but one of the weighting coefficients b s 
vanish. 

In the remainder of this section, we discuss some issues relevant to the experi- 
mental verification of these results. Consider a pump-probe measurement utilizing ei- 
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ther time-delayed photoionization or phase modulation with single-photon excitation, 
yielding a wave packet with p-state angular distribution. In time-delayed photoion- 
ization, the packet evolves after formation for a delay time t, whereupon it interacts 
with a probe pulse ionizing the atom. In the phase-modulation technique, the wave 
packet is excited by two identical laser pulses separated in time by a delay t. The 
interference of the two wave packets produces an ionization signal proportional to 
their overlap. 

Both these approaches permit measurement of an ionization signal displaying 



periodicities corresponding to the full/fractional revivals. Refs. fllOl , pTTJI describe the 
observation of a time-delayed photoionization signal for p-state radial wave packets 
in potassium for delay times up to approximately 800 psec. The measured ionization 
signal reveals both full and fractional revivals without any appreciable loss of the 



signal. Ref. p3| uses the phase-modulation technique for wave packets of rubidium 
with n ~ 46.5 and 53.3. Fractional revivals with periodicities as small as |T c i have 
been resolved. However, there was a decay of the signal envelope with increasing 
delay, and the signal was quite small after 300 psec. 

The limiting factor for the time delay in these experiments appears to be the abil- 
ity to maintain a good overlap of the two laser pulses in the interaction region as the 
delay is increased. For Rydberg states with n = 25 - 50, transitions induced by inter- 



actions with black-body radiation are not an issue in this context fl27| , p4| . Although 
the Rydberg-state lifetimes induced by black-body radiation are typically reduced by 
more than an order of magnitude compared to lifetimes induced by spontaneous emis- 
sion, they are still several orders of magnitude more than the superrevival time scale 
t sr . With a sufficiently low background pressure, transitions induced by collisional 
processes should not be a factor either. 



The time-delayed photoionization experiments in refs. [10, 11] maintained a good 
signal for time delays of up to 800 psec. This is long enough to allow detection of both 
full and fractional superrevivals. In the example of rubidium with n = 36, considered 
in Sec. VI, tf rac for q = 6 was found to be approximately 500 psec. Therefore, this 
full superrevival should be observable with existing experimental capabilities. The 
fractional superrevival at tf rac ~ 250 psec with T f * ac jt Tev should be measurable as 
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well. 

As n is increased, the number of autocorrelation peaks increases and the fractional 
superrevivals become more prominent. Values of n ~ 50 would require a delay line of 
3-4 nsec, and the overlap of the laser pulses would have to be maintained for this 
amount of time. Assuming this can be achieved, the observation of the superrevival in 
this case would show the return to the near classical state of a quantum wave packet 
after a time of over twelve hundred classical orbits. 
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APPENDIX A 



In this appendix, we show that at the times tf rac given in Eq. (10) the wave function 
ty(r,t) can be expanded as a sum of subsidiary waves of the form ip c \(f,t), with t 
shifted by a fraction of the period T c \. Section Al analyzes the phase periodicity of 
\l/(r, t) at the times if rac , while Sec. A2 obtains the phase periodicities of the subsidiary 
wave functions in the sum and matches them with the periodicity of \I/(r, t). We then 
examine the constraints on the times tf rac in A3 and the properties of the expansion 
coefficients in A4. 



1. Phase of Higher-Order Contributions in ^f(f,t) 

We first examine the periodicity at the times tf rac of the phase of the full wave 
function \I/(r, t). We are interested in the terms of higher order in k in Eq. (H) that 
result from the expansion of the energy in a Taylor series. 

Substituting the expression for t frac in Eq. (|10"D into \P(r, t) and neglecting terms 
in the exponential of order t rev /t ST < 1, we find that the second and third-order terms 
in the energy lead to a phase factor exp [27rz0fc], where 

k = ??*k?-™h>-*.k> . (32) 
Aq n q 
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We wish to find the minimum period I for shifts in k that leaves the phase exp [27ri0&] 
invariant. 

The analysis is simplified if we first impose a condition on m/n to reduce the 
general case to the same form as the one with n divisible by four. The definition for 
n in Eq. ( |TT| ) suggests we impose 

- = ^ (modi) , (33) 

which simplifies the expression for 9^ in the general case. The relatively prime integers 
m and n obeying Eq. (|33|) are then specified once the right-hand side of Eq. ( |33"D has 
been fully reduced. 

The choice (p3|) leads to a simplified form of the phase: 

6 k =**k>-*k> . (34) 
q q 

Define I as the minimum period in k that leaves 0& invariant (mod 1), i.e., I is the 
smallest integer satisfying O^+i = 6k (mod 1). Then, there are two possible solutions 
for I, both dependent on q: 

(q if q/9^0 (modi) , 

U/3 if q/9 = (mod 1) . 1 J 

2. Phase of the Subsidiary Waves ip c \(r,t) 

The second phase of interest results when the macroscopically distinct subsidiary 
waves i/) c i are shifted in time. We next evaluate the periodicity of this phase and 
match it to the one obtained above for \l/(r, t). 

Consider a shift At in time of the classical wave function ip c i(r,t) defined in Eq. 
(0). We take a shift of the form 

At = ^ rev + ^T cl , (36) 
b a 

where a and b are relatively prime integers, as are c and d. This shift generates a 
phase factor exp [— 2iri<j)k} in ^ c i(r,t), where 

2na c , ,„„, 
* t = — k + 2 k . (37) 
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To simplify this expression, we choose the relatively prime integers c and d to 
satisfy 

^ + 2 = (modi) , (38) 
where the first term is understood to be fully reduced. The phase 0& then reduces to 

- 8Va k . (39) 



3b 

To write the full wave function as an expansion in the subsidiary waves, the two 
sets of phases must have the same periodicity. We therefore impose the condition 
that the phase 4>k be periodic in k with the same minimum value I as the period of 
the phase 9k in *&(r, t). This requirement can be satisfied by an appropriate reduction 
of the fractional coefficient of k in Eq. ([39]). It implies that we can write 

8ria aN a 

iir = w = 7 • (40) 

The integer N is the product of all factors of 8rja that are also factors of I, and 
a = 8r]a/N is also an integer. Since a contains no factors that are also factors of I, a 
and I are relatively prime. With this reduction, <pk = ak/l. The minimum period for 
k in cj)k is therefore I, as desired. 

With the choices for c and d in Eq. fl38| ) and the reduction of 8r]a/3b in Eq. (|40|). 
we find that the time shifts At in ijj c i(f,t) become 

At = jT cl . (41) 



3. Constraints on tf rac 

The expansion of the wave function ty(r,t) in terms of the shifted wave functions 
ip c \ is given in Eq. (|13|) , and the b s coefficients are defined in Eq. (|14|). The time of 
formation of the subsidiary wave packets are the times tf rac given by Eq. (fLOf). The 
integers m and n must satisfy Eq. (p3|). We next prove that / and hence q must be a 
multiple of three. 

Examining the reduction of the ratio in Eq. (0), we see that the integer b satisfies 
the condition b = IN/ 3. Since all the factors of N are also factors of I, the only way 
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that b can be an integer is if I (and hence N as well) is a multiple of three. Therefore, 
the allowed values of I must all be multiples of three. Since I is given in Eq. fl3~5|) 
either as q for q/9 7^ (mod 1) or as q/3 for q/9 = (mod 1), it follows that q must 
also be a multiple of three. 

It turns out that the analysis of the wave-packet and autocorrelation-function 
periodicities given in Appendix B yields the further constraint p = 1. We therefore 
have the constraints on the times tf rac given by Eq. ( |TUD that q must be a multiple of 
three, p = 1, and m/n obeys Eq. (|33]). 



4. Properties of b s Coefficients 



The form of the complex coefficients b s multiplying the subsidiary wave functions 
ip c \ at the times tf rac , given in Eq. (0), reduce the expansion (13) to an identity. 
Substitute Eq. ([L4]) into the expansion (|13|) for \l/(r, t) at the times tf rac , and use the 
identity 



i-i 



7 E ex P 

6 s=0 \ 



-2ni 



.a(k-k') 



I 



(42) 



This yields the expression in Eq. (|8|) at times t ~ tf rac , up to terms in the exponential 
of order t Tev /t ST <^ 1. 

Note that the identity (|42|) holds only if a and I are relatively prime. This means 
that the reduction of the ratio of integers in Eq. ( [40]) is a necessary condition for the 
formation of the subsidiary wave packets at the time tf rac . Moreover, as discussed in 
A3, it leads to the restriction that q must be a multiple of three. 

If we multiply the expansion for b s in Eq. (|14D by its complex conjugate, sum over 
s, and use the identity in Eq. (f|2), we get 



l-i 

5>- 



i 



(43) 



This implies that the normalization of the wave function is maintained at the times 
^frac, as expected. 



APPENDIX B 

In this appendix, we examine the periodicities of the wave packet and the auto- 
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correlation function. Section Bl contains the proof that the wave packet \^(f, t)\ 2 is 
periodic at times t « if rac , with a period T frac that depends on i rev and T cl . Section 
B2 contains the corresponding analysis for the autocorrelation function. 



1. Periodicity of the Wave Packet 

The wave function \I/(r, t) through third order in the energy expansion is given at 
the times tf rac by Eq. fl8|). The issue of its periodicity can be addressed by considering 
the effect of shifting the time, 

e 

tfrac * ^frac ~jtvev , ('") 

where e and / are integers. 

We are interested in the smallest possible value of e/f that leads to periodicity in 
\I/(r, t). The appearance of periodic structures in \I/(r, t) upon time shifts of the form 
([44|) is possible only if the minimum periodicity of k in 6 k (cf. Eq. remains I and 
if the ensuing modified coefficients in the expansion ( |i3|) can be expressed in terms 
of the original coefficients b s . We have performed an analysis for arbitrary e and / 
along the lines of what follows in this subsection, from which we have shown that the 
minimum value of e/f is 

7-7, • 

For simplicity in what follows, we restrict our treatment here to this explicit value. 

Under the shift given by Eqs. ( fjlp and (15), the phase 6% of the higher-order 
time-dependent terms in \l/(r, t) becomes 

g, 3Q/P+1) k 2 _ P k S (46) 

q q 

where we have kept terms only up to order t rev /t ST <C 1. This phase continues to have 
the same minimum period I under changes in the summation index k as did the phase 
6k, i.e., 6' k+l = 6' k for the same values of / as in Eq. 

ui \n y i , tf r ac "T 

3, 



Next, consider the expansion of \l/(r, tf rac + ^t TCV ) as a sum of subsidiary wave 



functions tp c \ at the time tf rac + -t 



rev • 



3 , ti,, , _ ,2n sa s 



tf rac + -trev) = Y, K^cli?, tfrac + (— + ^)T cl ) . (47) 
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In this equation, we have used the replacement -t T 



2 n 



T cl , and the b' s are complex 



coefficients given by 



i-x 



^ = yE ex P 

' Ul._ 



fc'=0 



2m^k' 



exp {2ni9' k ,] 



(4* 



These coefficients are the same as those of Eq. (14) but with 6 k — >■ 6' k . 

We next show that the set of coefficients b' s in Eq. ( fi8| ) can be found in terms of 
the previous set b s , which are defined in Eq. First, we let k! — > k' — 1 in the 

definition of b' s . We get 



b' 





as" 




exp 


T 


exp 



2m 



3(r/p + 1) 



exp 



-2m- 



1-2 



y E ex P 



2m'^-k 



x exp 



2 / 3( W +l)-3p fc2 _ p^ 



exp 



2m 



6(r]p + 1) 3p N 



If we require that p = 1, then 

'3(^ + l)-3p 2 p 3 N 



exp 



2tt« 



g 



-A; - r -k' 
Q 



exp \2m6k 



(49) 



(50) 



where is defined in Eq. 



J 



We then define a quantity J as: 

3 if q/9 ^ (mod 1) , 
1 if q/9 = (mod 1) , 



(51) 



which permits us to write q = 3l/J. This gives 

.3(77 + 1) 



b' s = exp 



as 
T 



exp 



2m- 



x exp 



g 

2vri 



exp 



1 

-2m- 

g 



J-2 



- ^ exp [2m6 k ] 
1 k=-i 



'as J(2rj + 1) 
T + I 



(52) 



The quantity J(2i] +1)// can be rewritten as ax jl provided that an integer x 
exists that satisfies 

ax = J(2rj + 1) (mod /) . (53) 

It can be shown that an integer-valued equation of this form has a solution x if the 
greatest common divisor of a and I divides evenly into J(2r] + 1). Since a and I 
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are relatively prime, their greatest common divisor is 1, and this divides evenly into 
J(2?7 + 1). Thus, there always exists an integer x that satisfies Eq. (|3|), and we may 
replace (as + J(2r] + 1)) jl by a(s + x)/l. 

Since the exponential factors inside the sum fl52|) are periodic with period I, we 
can shift the summation over k back to the range k = to k = l — 1, We then obtain 



„ .as 
T 



K = exp 
Defining the phase factor 

exp [27ri6 s ] = exp 



exp 



2wi 



3(77 + 1) 



exp 



-2iri- 
Q 



.as 
I 



exp 



2n% 



3(77 + 1) 



exp 



-2m- 
Q 



we then have for p = 1 that the coefficients b' s are given by 



(54) 



(55) 



b' s = exp [27ri6 s ] b {s+x)r 



(56) 



where a; is a given integer. From this it follows that \b' s \ = |6( s+x .) mod J. We use this 
relation to prove that the motion of the wave packet is periodic. 

The probability density for the wave packet at the time tf rac + H Tev is the absolute 



square of the expression fl47|). Since the subsidiary wave packets are by construction 
localized in space and distributed evenly along the orbit, they have little or no overlap. 
This permits us to ignore the cross terms in the absolute square of the wave function. 
Substituting the expression fl56|) for b' s into Eq. ( pET|) and taking the absolute square 
gives 



3 2 

ty{7'. tfxac + — ^rev)| 



,2n sa. 



Yl lfys+x) modi lVcl(^frac + (— + -r)T< 

q i 



cl 



(57) 



1 s=0 

Since the integer x always exists, we can set s' — (s + x) mo d/ and reorder the terms 
to get 



z-i 



rev j 

q 



J2 l^'HV'cl^tfrac + (— - ^r)T d + ^-T cl ) 

5Tn q I I 



(5f 



Introduce the irreducible ratio u/v of integers u and v satisfying 

u 2n ax 

— = (mod 1) . 

v q I 



(59) 
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Using Eq. (|TTJ) of Sec. IIIA and Eqs. fl5T|) and fl53|), this definition of u/v may be 
rewritten as 

u 2(r? + A) -3 



v q 
Equation fl58|) then becomes 



(mod 1) . (60) 



3 u 
|*(r,tfrac + -t ICV )\ 2 = |*(r,tfrac + -^ci)| 2 , (61) 

where we have used Eq. ( |i~3| ) and the assumption that the subsidiary wave packets 
do not overlap. 

Since T c \ <C t rev , the expansion in Eq. ( jl3|) holds for shifts in time of order T c \. 
We may therefore subtract ^T cl from t on both sides of Eq. ( |6"ID to obtain 

3 it 

l*^, tfrac + "trev ~ -T d )| 2 = | (f, £ frac ) | 2 . (62) 



Finally, we can define the period Tf rac as in Eq. fliq) , where the integers u and t> 
satisfy Eq. (|60|) . 

We have thus proved that at the times t ~ tf rac , the wave packet is periodic with 
the period Tf rac . 



2. Periodicity of the Autocorrelation Function 

A similar method can be used to show that at the times tf rac the absolute square 
of the autocorrelation function is periodic with the same period Tf rac . The autocorre- 
lation function at time tf rac is A(tf Tac ) = ($(r, 0)|\l/(r, tf ra c))- The initial wave function 
\l/(r, 0) is ipd(r, 0). Shifting tf rac to tf rac + H TCV and using Eqs. (f47|) and (j56|) and the 
definition of u/v in Eq. flBUD yields 

3 i_1 u as' 

A(t irac + -t IeY ) = E exp [27iiQ s >] b s > (V>ci(0)|^ c i (tfrac + -T c \ + —j-T c \)) , (63) 
q s / =0 v l 

where we have suppressed the r-dependence. The phase exp [27riO s /] is given in Sec. 
Bl. Note that the substitution s' = (s + x) mo di has been made. 
Taking the absolute square of A(tf rac + ^t Tev ), gives 

3 l ~ 1 u as' 

\A{t h&c + -t rcv )| 2 = M 2 K^cl(0)|^cl(tfrac + -T d + ^T cl ))| 2 

9 T^o v I 
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I— 1 U OlS 

+EE ,± -n eX P M©*' - ©•)] ^ (^cl(tfrac + -T d + - r T c l)|V'cl(0)) 

X (^cl(O)l^cl(tfrac + ~T d + — T cl )) . (64) 
V I 

The first term on the right-hand side is just | A(tf rac + |T c i)| 2 . Since s ^ s' in the 
second term and since we assume that the subsidiary wave functions do not overlap 
appreciably, one of the inner products in the double sum must vanish. We can 
therefore neglect this term. 

Subtracting ^T c \ from both sides produces 

3 u 

|^frac+-trcv--T c l)| 2 = |A(t frac )| 2 . (65) 

q v 

This proves that the square of the autocorrelation function is periodic at the times 
t frac with period T frac given by Eq. (jTJ). 

APPENDIX C 

In this appendix, we prove that at the times tf rac the SQDT wave packet \l/(r, t) for 
an alkali- metal atom can be written as an expansion over subsidiary packets ip c \{r, t). 
In analogy with the method for hydrogen, we proceed by matching the period of the 
phase of the higher-order contributions in \l/(r, t) with the period of the phase induced 
by shifting ip c \(r, t) by a fraction of its period T c *. We then determine the coefficients 
in the expansion and the constraints on the times tf rac . Lastly, we show that the wave 
packet and autocorrelation function are periodic with a period T f * ac and we determine 
the allowed values of T f * ac . 

Consider ^(f,t) at the times t£ rac . An additional phase exp [27rz0j£] due to the 
higher-order corrections appears relative to the linear term, where 

ei = ^e-™e- p -e . (ee) 

4q n q 

As usual, terms of order t* ev /t* r <C 1 in the exponential have been neglected. 

To simplify the analysis as much as possible, we choose the integers m and n so 
that 91 reduces to the value of 6 k defined for hydrogen in Eq. (|34T). This imposes 

m 3p{Xu-fi) 

— = -: (mod 1) , (67) 

n Aqu 
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where the right-hand side is understood to be fully reduced. With this choice, the 
phase 9 k equals the phase 9 k for hydrogen given in Eq. fl3"3] ) . It therefore has the same 
minimum period I under shifts in the summation index k as for hydrogen, given in 
Eq. ©. 

Note that if fi = 0, the definition (B^) for m/n reduces to the form in Eq. (|33| ) for 
hydrogen. The additional shift for // 7^ comes from the fractional part of n* due to 
the quantum defect, assuming the laser is on resonance. We therefore see that the 
allowed times tf rac are not simply obtained by the scaling transformations t sr — > t* r 
and t rev — > t* ev in the definition of £f rac . In addition to the scaling, the fraction m/n 
is shifted by an amount that depends on the quantum defect. 

We next consider the phase 4>* k induced by shifting the classical waves ip c i(r,t) by 
an amount 

At = ^ev + ^ , (68) 

where a, b and c, d are pairs of relatively prime integers. The phase generated is 

2n*a, c, , , 

« = — k + -k . (69) 

Substituting the definition of fi* given in Eq. (|28|) , we choose the integers c and d 
to satisfy 

2(\v — a) a c „ , , , 

3^ + ^ ( } ' ( } 

With this condition, (p* k reduces to the expression ( ff^D for the phase 0^ in hydrogen. 

Following the definitions in Eq. ( ^0|) , cf) k can be written as a /Ik, where a and I are 
relatively prime integers as defined in Sec. IIIC The time shifts At are then given by 

At = j23 • (71) 

The minimum period of <p* k under shifts in the summation index k is again /, which 
matches the period of 9 k . 

Since the phases induced by the shifts At in ip c i(r,t) have the same period as 
the higher-order contributions to the time-dependent phase in \l/(r, t) at the times 
tf rac , we may use the set 4) c \{r, t* h&c + ^fT* { ) with s = 0, 1, — 1 as a basis for an 
expansion of \&(F, tf rac ), and we obtain Eq. fl3"0p. Since the phases 6*^ and <p* k have the 
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same formal structure as those in hydrogen, the coefficients b s may be taken to have 
the same form as in Eq. flU]). The proof that this expansion is valid follows as before, 
using Eq. (g^). 

The allowed values of q in the definition (p7|) of tf rac are restricted to multiples 



of three. The proof follows that in hydrogen. Since the phase (j)* k reduces to <pk for 
hydrogen, and since Eq. ( pEUp holds, the quantity IN/3 must again be an integer. This 
is true only if I and hence q are multiples of three. 

We next prove that the wave packet and absolute square of the autocorrelation 
function are periodic for times near if rac and we determine the periodicity T f * ac . The 
procedure is similar to the corresponding proofs for hydrogen. 

Consider the wave function \P(r, £f rac ) under a shift in time t* h&c — > t* h&c + jt* ev . It 
can be shown that the minimum shift that leads to periodic behavior is e/f = 3/q. 
The wave function at this shifted time can be written as an expansion in ip c \ with 
coefficients b' s given in Eq. (¥3). In parallel with the hydrogenic case, for p = 1 we 
find the relation \b' s \ = \bt s + x ) mod i\ f° r some integer x that can be shown to exist. 

Taking the squared modulus of the wave function and using the assumption that 
the subsidiary waves do not overlap, we deduce 

l*(^^ a c + ^ac)| 2 =l^(^^ac)| 2 , (72) 

with the period T f * ac = ^t* cv — ^T c *. The integers u and v satisfy a relation that 
depends on the fractional part of the quantum defect: 

„ = 2(„ + A)-3_2, 
v q qv 

Equation ( ff2|) shows that near the times tf rac the wave packet is periodic with period 

rri* 
frac " 

Using a similar technique, we have proved that the square of the autocorrelation 
function is also periodic at the times £f rac , with the same period T f * ac . 
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FIGURE CAPTIONS 



Fig. 1: The absolute square of the autocorrelation function for a hydrogenic Rydberg 
wave packet with n = 319 and a = 2.5 is plotted as a function of time in 
microseconds. The plots span times of order of the revival time t rcv ~ 1.05 
//sec. (a) t — to £ — 0.6 //sec, (b) t = 0.6 to £ — 1.2 //sec. 

Fig. 2: The absolute square of the autocorrelation function for a hydrogenic Rydberg 
wave packet with n = 319 and a = 2.5 is plotted as a function of time in 
microseconds, (a) t — to t — 15 //sec, (b) £ = 15 to t = 30 //sec, (c) £ = 30 to 
£ = 45 //sec. 

Fig. 3: Unnormalized circular wave packets for hydrogen with n = 319 and a = 
2.5. Cross-sectional slices of the wave packet in the plane of the orbit and for 
r — (r) — \n(2n + 1) are plotted as a function of the azimuthal angle <f> in 
radians, (a) £ = 0, (b) £ = £ rcv , (c) £ = |£ sr - §£ rcv , (d) £ = |£ sr + |£ rcv - \T ch 

(e) £ 9^ sr 4^rcvj (f) £ 9^ sr I2^ rev 9"^°'' (§) ^ 12 ^ sr X6^ rcv ' ( ' 

t 12 ^ sr ig^rev ^2-^~ch W t 15 ^ sr 20^ rcv ' 0) ^ l~5^ sr 20^ rcv Y5 -^cl ? 0^) 

i — J-/ — 1/ (V\ i — -L+ -4- -L+ _ II T 

18 sr 8 l rev) vv ^ 18 sr 24 rev 18 c '' 

Fig. 4: Unnormalized radial wave packets for hydrogen with n = 36 and o = 1.5. The 
radial probability density is plotted as a function of r in atomic units, (a) £ = 0, 

(b) £ i^Tcl) (c) £ £rcv "I - 2"^'' (^) ^ 6^ sr 6"^°'' ^ 6^ sr 2^ rcv 3~^e'' 

Fig. 5: The absolute square of the autocorrelation function for a Rydberg wave packet 
for rubidium with n = 36 and a — 1.5 is plotted as a function of time in 
nanoseconds. The p-state quantum defect is 5(1) = 2.65. 

Fig. 6: Unnormalized radial wave packets for rubidium with n = 36 and a = 1.5. The 
p-state quantum defect is 5(1) = 2.65. The radial probability density is plotted 
as a function of r in atomic units, (a) £ = 0, (b) £ = \T* h (c) £ = £* ev + j^T* h 

(J\ j- l-L* 27 j.* I 2rp* i \ . ]..* , 53 J-* _ 13 rp* 

V / 6 sr 160 rev ~r 5 J cl> l e J 1 ~~ 6 sr ~r 160 rev 60 cl" 
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